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1.0  INTRODUCTION 


Atmospheric  turbulence  affects  the  performance  of  missile  defense  platforms,  such  as  the 
Airborne  Laser  (ABL).  To  understand  the  results  of  laser  propagation  tests  on  both 
horizontal  and  inclined  propagation  paths,  and  to  validate  atmospheric  turbulence 
prediction  models,  the  vertical  profile  of  the  refractive  index  structure  characteristic 
Cl{h)  is  required.  The  astronomical  and  adaptive  optics  communities  have  established 

00 

several  optical  techniques*  for  measuring  path-integrated  values  of  =  ^Cl{h)dh  or  the 

0 

Fried  parameter,  ,  but  these  methods  do  not  provide  path-resolved  information.  In 
order  to  remotely  sense  a  vertical  profile  of  ,  a  true  lidar  is  required. 

A  number  of  lidar  techniques  for  measuring  C^(A)  have  been  proposed.  These  methods 
were  based  on  turbulent  broadening  of  a  laser  beam,^"^  the  residual  turbulent  scintillation 
effect,^  and  laser  beam  motion.®  Measurements  of  atmospheric  Cl  values  with  a  bistatic 
CW  lidar  were  reported  in  Ref.  3.  However,  despite  these  developments,  there  is 
currently  no  accepted  lidar  method  for  measuring  profiles  of  Cl  {h) . 

The  differential  image  motion  (DIM)  method  has  already  been  established  in  astronomy 
as  the  best  way  to  measure  the  Fried  parameter,  ,  but  all  of  the  previous  work  used  a 
star  or  a  laser  at  one  end  of  an  optical  path  and  a  receiver  at  the  other.  For  this  reason, 
only  path-integrated  quantities  such  as  can  be  measxured  with  this  method. 

The  Applied  Technologies  Inc.  and  Georgia  Tech  Research  Institute  team  proposed  a 
new  concept,  DIM  lidar,  for  measuring  vertical  profiles  of  Cl{K) .  A  lidar  has  the 

transmitter  and  receiver  at  the  same  end  of  the  optical  path,  and  it  uses  a  pulsed  laser, 
atmospheric  backscatter,  and  a  range-gated  receiver  to  probe  the  atmosphere  at  various 
ranges.  Phase  I  of  the  SBIR  was  an  investigation  of  this  concept. 
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2.0  PHASE  I  TECHNICAL  OBJECTIVES  AND  MILESTONES 


The  overall  objective  of  Phase  I  was  to  demonstrate  the  technical  feasibility  of  the  DIM 
lidar  concept  by  developing  and  testing  a  hard-target  analog  of  this  lidar. 

The  Phase  I  work  plan  included  six  tasks: 

•  Conceptual  design; 

•  Performance  analysis; 

•  Hardware  and  software  development; 

•  Field  test; 

•  Data  analysis;  and 

•  Scientific  paper  and  final  report. 

The  work  was  planned  for  nine  months.  A  mid-term  review  and  a  briefing  were 
scheduled  on  September  2,  1999,  and  delivery  of  the  final  report  was  scheduled  in 
February,  2000. 

3.0  ACCOMPLISHMENTS 

All  planned  tasks  were  accomplished.  The  results  from  Phase  I  were  presented  at  a 
briefing  at  the  sponsor  facility  on  September  2,  1999.  Hardcopies  of  the  briefing 
materials  were  delivered  to  the  Technical  Monitor.  Additional  work,  which  includes  the 
analysis  of  the  effects  of  a  speckle  phenomenon  and  scintillation  on  the  performance  of 
the  DIM  lidar,  was  documented  in  a  Technical  Memorandum  dated  September  1999, 
which  was  also  delivered.  The  accomplishments  for  each  task  are  described  in  the 
following  sections. 
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3.1  Conceptual  Design 


3.1.1  DIM  lidar  concept 

The  concept  of  the  DIM  lidar  is  illustrated  in  Fig.l.  This  lidar  employs  common 
transmit/receive  optics  and  a  range-gated  focal  plane  array.  The  lidar  operates  as  follows: 
pulsed  laser  beam  is  transmitted  through  a  common  transmit/receive  aperture  and  is 
focused  at  the  desired  measurement  distance.  Scattered  laser  irradiance  is  received 
through  two  spatially  separated  sub-apertures.  The  motion  of  two  images  of  the  laser 
beam  in  the  scattering  volume  is  measured  in  the  image  plane.  The  statistics  of  the 
differential  motion  of  the  two  images  relates  to  the  strength  of  optical  turbulence  along 
the  path,  which  is  characterized  by  the  refractive  index  stracture  characteristic  .  By 
measuring  the  differential  motion  of  a  reflected  laser  spot  over  various  distances,  a  profile 
of  is  estimated. 
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3.1.2  Theoretical  underpinnings  of  the  DIM  lidar 


An  analytical  expression  for  the  differential  image  motion  variance  is  known  for  the 
conventional  DIM  technique.^  For  the  lidar  configuration  it  was  established  under  this 
task.  In  the  lidar  scheme,  an  optical  wave  propagates  through  the  turbulent  medium 
twice:  fi-om  the  lidar  to  the  scattering  volume,  and  back  to  the  receiver.  For  a  common 
transmit/receive  configuration,  the  image  displacements  can  be  presented  in  the  form 


<PX=(PB-(PX,R,  (1) 

9i=(Pb-^i,r,  (2) 


where  and  (p^  are  the  eentroid  positions  of  the  two  images,  is  the  displacement  due 
to  the  wavefront  tilt  of  the  transmitted  beam,  and  <p^^  and  are  the  displacements  due 

to  the  wavefront  tilts  of  the  reflected  wave  in  the  two  receiving  apertures  The  minus 
signs  in  equations  (1)  and  (2)  take  into  account  that  the  transmitted  beam  and  the 
reflected  wave  propagate  in  opposite  directions.  By  subtracting  Eq.  (1)  from  Eq.  (2)  one 
can  show  that  the  differential  image  motions  on  one-way  and  two-way  propagation  paths 
coincide  with  each  other,  i.  e.  (p^  =■  <Px—<P2-  (Pi^R  ~  (Px^r  »  and  so  the  differential  image 
motion  variance  for  a  lidar  is  given  by 

O-rf  =  {9])  =  {(92.R  -  <Px,Rf)-  (3) 

The  implication  is  that  one  can  use  the  same  expression  derived  for  the  eonventional  one¬ 
way  scheme’  for  the  lidar  configuration.  This  expression  has  the  form 


11.6[l-0.694(i/ ///)''"]’ 


(4) 


_  L 

where  =  Jc’  (x)(l -x!  Lf'^dx  ,L\s\hQ  distance  between  the  lidar  and  the  scattering 
0 
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volume,  /  is  the  focal  length  of  the  receiving  telescope,  d  is  the  sub-aperture  diameter, 
and  //  is  the  separation  between  the  sub-apertures.  Eq.(4)  takes  into  account  the 
spherical  divergence  of  the  reflected  wave.  For  a  horizontal  path  where  C]  is  constant, 
cl  =  {3/8)LC^  .  Eqs.  (3)  and  (4)  is  also  valid  for  a  bistatic  lidar  configuration. 

The  absolute  value  of  the  image  motion  variance  in  each  receiving  channel  of  the  DIM 
lidar  for  common  transmit/receive  optics  is  given  by* 

0 

where  is  the  transmitting  aperture  diameter,  and  >  d .  For  a  bistatic  configuration, 

where  the  wavefi'ont  tilt  of  the  transmitted  beam  and  the  reflected  wave  are  uncorrelated, 
this  equation  has  the  form 


(«)5  )  =  4.04(a-''>  +d-''>)]cl  «)(1  -ILy'‘d4.  (6) 

0 

According  to  Eq.  (5),  when  the  transmitter  diameter  is  equal  to  the  receiver 
diameter  a^-^d,  the  centroid  variance  equals  zero^^j^2^ 0.  This  is  caused  by 

reciprocity,  i.e.  the  fact  that  the  transmitted  beam  and  the  reflected  wave  propagate 
through  the  same  turbulent  eddies.  However,  even  for  common  transmit/receive  optics 

the  absolute  variances  (^1^2)  ^  non-zero  value  when  a^>d .  Besides,  fi'om  Eqs. 

(l)-(3)  it  is  seen  that  since  the  differential  image  motion,  q>^ ,  does  not  depend  on  the 
beam  wander  component,  ,  so  reciprocity  does  not  affect  the  performance  of  the  DIM 
lidar. 

The  DIM  method  has  numerous  advantages  over  other  proposed  techniques.  Because  it 
measures  the  differential  motion  of  images  and  not  their  sizes,  it  is  fundamentally 
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immune  to  any  effects  that  cause  image  blurring,  such  as  vibration,  laser  instability, 
diffraction  on  the  transmitting  and  receiving  apertures,  environmental  changes  in  the 
system  MTF,  and  aerosol  image  blurring.  It  is  not  affected  by  reciprocity,  and  the  DIM 
technique  has  been  shown  to  be  immune  to  variations  of  the  outer  scale  of  turbulence. 
This  effect  is  considered  in  the  next  section.  For  all  these  reasons,  the  DIM  technique  is 
much  more  robust  than  methods  based  on  laser  beam  broadening  and  beam  motion. 

In  addition,  the  proposed  lidar  is  based  on  direct  detection  (rather  than  coherent 
detection),  which  means  that  it  can  be  built  largely  from  commercial,  off-the-shelf 
components,  and  the  optical  assembly  will  not  have  a  requirement  for  extreme  stability. 
The  astronomical  and  adaptive  optics  communities  have  already  established  the 
technology  base  for  the  DIM  lidar.  Development  of  a  prototype  would  be  low-risk  and 
cost-effective. 

3.1.3  Effect  of  the  outer  scale  of  turbulence 

An  analysis  of  the  effect  of  the  outer  scale  of  turbulence  on  the  DIM  technique  is 
presented  in  this  section.  In  an  experiment  performed  at  the  Starfire  Optical  Range  (SOR) 
at  Kirtland  Air  Force  Base’  three  optical  methods  including  a)  the  DIM  technique,  b)  the 
long-exposure  PSF,  and  c)  the  short-exposure  PSF  methods  were  compared.  The  Fried 
parameter  measured  with  these  three  sensors  is  shown  in  Fig.  2.  The  data  from  the  DIM 
and  long-exposure  PSF  sensors  are  marked  as  ATMOS  and  MTF,  respectively,  whereas 
the  data  from  the  short-exposure  PSF  sensor  are  marked  as  SHORT  EXP. 
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2B  Oct  1996,  Starfirc  Optical  Range 
r^,  zenith  corrected,  500  nm 


Figure  2.  The  Effect  of  Outer  Scale 

The  comparison  revealed  that,  for  some  conditions  significant  differences  (up  to  a  factor 
of  2)  between  the  readings  from  the  three  sensors  has  been  observed.  An  analysis 
performed  by  researchers  at  Georgia  Tech  showed  that  the  three  methods  have  different 
sensitivities  to  the  variations  of  the  outer  scale  of  turbulence,  Lo,  and  the  DIM  method  is 
the  least  sensitive  to  this  effect.  This  is  because  large-scale  wavefront  aberrations  caused 
by  the  outer  scale  of  turbulence  are  cancelled  in  the  differential  motion.  This  analysis 
was  performed  by  using  the  von  Karmann  power  spectrum  of  turbulence,  which  includes 
a  finite  outer  scale,  Lq  . 

3.2  Performance  Analysis 
3.2.1  3.2.1.  Turbulent  calculations 

Under  this  task  a  performance  analysis  of  the  DIM  lidar  was  performed,  including 
quantitative  characterization  of  the  laser  spot  size  in  the  image  plane  and  the  amount  of 
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motion  of  this  image  spot.  This  analysis  also  includes  characterization  of  the  signal-to 
noise  ratio,  and  the  centroid  error  caused  by  a  finite  number  of  photons.  Parameters  of  the 
laser  beam  and  reflected  spot  in  the  image  plane  were  estimated  for  four  propagation 
scenarios:  1)  a  proof-of-concept  test  over  300  m,  2)  tests  at  the  ABL  advanced  Concept 
Test  Bed  (ACT),  3)  profiling  of  at  astronomical  sites,  and  4)  measurements  of 

aboard  the  ABL  platform.  Since  the  above  scenarios  include  different  measurement 
ranges,  which  vary  from  300  m  up  to  300  km,  different  laser  sources  were  considered  for 
each  case.  A  He-Ne  25  mW  laser  was  considered  for  the  proof-of-concept  test  over  300 
m.  A  Q-switched,  doubled  Nd:YAG,  operating  at  532  nm  and  transmitting  10  W  laser 
power  was  planned  for  the  ABL  ACT  test.  An  eye-safe  UV  laser  operating  at  351nm  was 
considered  for  the  commercial  astronomical  application.  Finally,  a  pulsed  laser  operating 
at  1.06  pm  with  an  energy  1  J  per  pulse  and  a  pulse  repetition  rate  of  1  kHz  was  planned 
for  the  ABL  scenario. 

The  estimated  parameters  include  the  coherence  radius,  or  Fried  parameter,  r^,  long- 
exposure  beam  size,  short-exposure  beam  size,  (p^j.  (which  determines  an 

instantaneous  spot  size  in  the  image  plane),  beam  wander  standard  deviation,  cr^,  and 
rms  differential  motion,  cr^  .  The  long-exposure  angular  size  of  the  beam  is  given  by  the 
known  formula^ 

where  L  is  the  distance,  is  the  initial  radius  of  curvature  of  the  beam,  is  the  radius 
of  the  transmitting  aperture,  k  =  lnlX  is  the  wavenumber,  and 

L 

P(~  =  (1.45^^  |c^(z)(l  -z/ Ly^dzY^'^  is  the  coherence  radius  of  a  spherical  wave.  For 

0 

cl  =  const,  this  radius  relates  to  the  coherence  radius  of  a  plane  wave,  p^ ,  by 
P(,  =  l.Spo  and  to  the  Fried  parameter,  Tq,  by  =  l.lp,- . 
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The  short-exposure  radius  of  the  beam,  according  to  Ref.  10,  relates  to  the  long- 
exposure  beam  radius,  ,  and  rms  beam  wander,  ^  ,  by  the  equation 

~^SE 

Consequently,  the  short-exposure  beam  radius  is  The  two-axis 

beam  wander  variance  for  a  focused  beam,  has  the  form'® 

Table  1  presents  the  estimated  parameters  of  the  transmitted  beam  and  reflected  laser  spot 
in  the  image  plane  for  the  proof-of-concept  test.  The  expected  values  of  Cl  are  from 
10’'‘’w”^^^to  the  transmitter  diameter  is  la^  =  0.09m ,  the  sub-aperture 

diameter  and  the  separation  are  d  =  0.\\m  and  //  =  0.1  Iw ,  respectively.  It  is  seen  that 
the  short-exposure  beam  size  is  in  the  range  from  3.5  /jrad  to  39  jM'ad ,  whereas  rms 
differential  motion  is  in  the  range  from  2.9  /jrad  to  20  ftrad .  Importantly,  that  rms 
differential  motion  is  comparable  to  the  short  exposure  size  of  the  beam,  and 
consequently,  to  the  instantaneous  spot  size  in  the  image  plane.  This  suggests  that  the 
differential  image  motion  can  be  accurately  measured.  The  proof-of-concept  test 
performed  under  task  4  confirmed  this  conclusion. 

Tables  2,  3,  and  4  present  estimated  parameters  of  the  transmitted  beam  and  reflected 
laser  spot  in  the  image  plane  for  the  ABL  ACT  scenario.  The  selected  values  of  Cl  ( 

2-10"‘’w'^^%  6-10“”m"^'^and  2-10"‘®w"^'^),  correspond  to  the  values  of  the 
Fried  parameter  of  =20  cm,  10  cm,  and  4  cm  measured  at  the  wavelength  of  1 .06  jum  on 

the  50  km  propagation  path.  The  absolute  values  of  the  estimated  rms  differential  motion 
and  short-exposure  beam  size  in  this  case  are  smaller  than  that  for  the  proof-of-concept 
test,  but  they  are  comparable  to  each  other.  Therefore,  the  image  motion  of  a  reflected 
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spot  can  be  accurately  measured  for  the  ABL  ACT  test  as  well.  The  estimates  given  in 
Tables  2,  3,  and  4  would  be  useful  for  the  design  of  the  DIM  lidar  for  the  ABL  ACT  test 
site. 

Table  5  shows  system  parameter  estimates  for  the  astronomical  scenario.  An  eye-safe  UV 
laser  was  considered  in  this  case.  The  path-integrated  values  of  C/ ,  which  correspond  to 
values  of  the  Fried  parameter  in  the  range  from  5  cm  to  20  cm,  were  selected.  These 
values  of  correspond  to  a  500  nm  wavelength  and  are  corrected  to  the  zenith. 

Finally,  tables  6  and  7  present  system  parameter  estimates  for  the  ABL  scenario.  In  this 
case  a  1  m  telescope  and  sub-apertures  of  0.3  m  in  diameter  with  0.3  m  separation  were 
considered.  The  propagation  ranges  were  100  km  and  300  km.  The  estimates  in  Tables 
5-7  also  show  that  the  image  motion  can  be  accurately  measured. 


TABLE!.  PROOF-OF-CONCEPT  TEST 

L  =  300M,  X=  0.63xl0’*M,  2a„  =  0.09M,  d  =  O.llM,  n  =  O.llM 


Structure 

Parameter 

C„',  M-“ 

10'^ 

5x10*'’ 

10'^ 

5x10-'^ 

Coherence 

Parameter 

To, cm 

5.5 

2.1 

1.4 

0.56 

Long-Exposure 

Beam  Size 

(peff,liRad 

5 

12 

17.6 

46.6 

RMS  Beam 

Wander 

Qc,  ,pRad 

3.6 

8 

11.4 

25.6 

Short-Exposure 

Beam  Size 

(pST,tiRad 

3.5 

4.5 

13.4 

39 

RMS  Differential 
Motion 

CTd,|xRad 

2.9 

6.5 

9.1 

20.4 
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TABLE  2.  AIRBORNE  LASER  ADVANCED  CONCEPTS  TESTBED 


Focused  beam 


=  1) ,  0.532xl0  *M,  2ao  =  0.2M,  d  =  0,16M,  =  0.16M,  =  2xlO  *’]Vr^'^ 


Range 

L,  km 

10 

20 

30 

40 

50 

Coherence 

Parameter 

To, cm 

22 

15 

12 

9.9 

8.6 

Long-Exposure 

Beam  Size 

(peff,|iRad 

1.2 

1.6 

1.9 

2.6 

2.4 

RMS  Beam 

Wander 

Gc,  ,|LiRad 

0.8 

1.1 

1.4 

1.6 

1.8 

Short-Exposure 

Beam  Size 

1.1 

1.3 

2.1 

1.6 

RMS  Differential 
Motion 

ad,|iRad 

0.7 

1.0 

1.2 

1.4 

1.6 

TABLE  3.  AIRBORNE  LASER  ADVANCED  CONCEPTS  TESTBED 


Focused  beam  (—  =  1) ,  0.532x1 0  'TVI,  2a„  =  0.2M,  d  =  0.16M,  n  =  0.16M, 

Fi 


Range 

L,  km 

10 

20 

30 

40 

50 

Coherence 

Parameter 

To, cm 

12 

7.7 

6 

5.1 

4,4 

Long-Exposure 

Beam  Size 

(PeflvliRad 

1.9 

in 

3.3 

3.9 

4.5 

RMS  Beam 

Wander 

Gc,  ,^Rad 

1,4 

2.0 

2.4 

2.8 

3.2 

Short-Exposure 

Beam  Size 

<PsT,M-Rad 

1.25 

1.8 

2.3 

2.8 

3.2 

RMS  Differential 
Motion 

CTd,nRad 
_ 1 

1.21 

1.7 

2.1 

2.4 

2.7 
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TABLE  4.  AIRBORNE  LASER  ADVANCED  CONCEPTS  TESTBED 


Focused  beam  (—  =  1) ,  0.532x10  ‘M,  lap  =  0.2M,  d  =  0.16M,  ii  =  0.16M,  C„^  =  2x10  ‘®M  “ 

Fl 


Range 

L,  km 

10 

20 

30 

40 

50 

Coherence 

Parameter 

To, cm 

5.7 

3.7 

3.0 

2.5 

2.2 

Long-Exposure 

Beam  Size 

(Peff.nRad 

3.5 

5.3 

6.7 

7.6 

9.1 

RMS  Beam 

Wander 

CTc,  ,^Rad 

2.6 

3.6 

4.4 

5.1 

5.7 

Short-Exposure 

Beam  Size 

(psT,nRad 

2.5 

3.9 

5.0 

5.7 

7.1 

RMS  Differential 
Motion 

Od,nRad 

2.2 

3.1 

3.8 

5.0 

TABLES.  ASTRONOMY 

L  =  20km,  0.351x10-*]^,  2ao  =  0.5M,  e  =  90%d  =  0.16M,  |i  =  0.1 6M 


Structure 

Parameter 

c 

n 

2.2x10-'^ 

7x10-'^ 

2.2xl0-'^ 

10" 

Coherence 

Parameter 

ro,cm 

13.1 

6.5 

3.3 

1.3 

Long-Exposure 

Beam  Size 

(Peff,^Rad 

1.0 

2.1 

4.1 

10.1 

RMS  Beam 

Wander 

Op,  ,^Rad 

0.8 

1.4 

2.4 

5.2 

Short-Exposure 

Beam  Size 

(psT,liRad 

0.7 

1.6 

3.3 

8.7 

RMS  Differential 

Motion 

ad,nRad 

0.7 

1.3 

2.3 

5.0 
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TABLE  6.  AIRBORNE  LASER  SCENARIO 


L  =  lOOkm,  X=  1.06xl0  *M,  2ao  =  IM,  d  =  0.3M,  ^  =  0.3M 


Structure 

Parameter 

10-'* 

5x10-'* 

10-’’ 

5x1 0-'^ 

Coherence 

Parameter 

To, cm 

79.1 

30.2 

20.0 

7.6 

Long-Exposure 

Beam  Size 

(Peff,nRad 

0.6 

1.4 

2.1 

5.4 

RMS  Beam 

Wander 

ac,  ,|iRad 

0,46 

1.0 

1.45 

3.25 

Short-Exposure 

Beam  Size 

cpsT.liRad 

0.4 

0.95 

1.5 

4.3 

RMS  Differential 

Motion 

ad,nRad 

0.4 

0.92 

1.3 

2.9 

TABLE?.  AIRBORNE  LASER  SCENARIO 

L  =  300kin,  X=  1.06xl0  ®M,  2ao  =  IM,  d  =  03M,  =  0.3M 
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3.2.2  Lidar  calculations 


The  uncertainty  in  the  measurement  of  an  image  centroid  position  due  to  photon  statistics 
is  given  by 


(10) 


where  SJs  the  position  uncertainty,  cr^ is  the  image  width,  Nj^is  the  number  of  laser 

photons  in  the  image,  and  A^^is  the  number  of  sky  background  photons  in  the  image. 

The  number  of  laser  photons  was  calculated  by  using  the  lidar  equation,  which  can  be 
expressed  as 


where  N  is  the  number  of  photons  received. 

No  is  the  number  of  photons  transmitted  in  each  laser  pulse, 

A  is  the  receiver  area  R  is  the  range  (m), 

k  is  the  optical  efficiency  (dimensionless), 
c  is  the  speed  of  light  ^  x  lO^w/5^, 

r  is  the  sampling  interval  (s), 
p  is  the  backscatter  coefficient  (m  ^sr  ,  and 

a{r)  is  the  atmospheric  extinction  coefficient  (m 


(11) 


The  number  of  photons  transmitted  is  related  to  the  energy  per  pulse  E  in  Joules  by 

No  =  AE/hc,  (12) 

where  X  is  the  wavelength  (in  meters),  h  is  Planck’s  constant,  6.63  x  10“^'*  Js ,  and  c  is 
the  speed  of  light. 
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The  number  of  sky  backgroimd  photons  was  calculated  from  the  equation 


P=LAQB,  (13) 

where  P  is  the  received  power  (W), 

L  is  the  spectral  radiance  of  the  sky  (w  Im^  -  sr-  micron^, 

A  is  the  receiver  area 

Q  is  the  receiver  solid  angle  (sr),  and 
B  is  the  receiver’s  optical  bandpass  (microns). 

The  solid  angle  Q(5r)  can  be  calculated  from  field  of  view  angle  ^(rad)  by  using  the 
relation 


Q  =  (14) 

The  number  of  sky  background  photons  that  will  be  received  by  the  LIDAR  during  the 
sampling  interval  r  is  given  by 


Ng  =  PrAlhc  ,  (15) 

so  by  substituting  the  expression  in  Eq.  (13)  for  P  we  find 

=  LAQBrX/hc .  (16) 

In  order  to  estimate  centroid  position  uncertainties,  the  numbers  of  photons  and 
were  calculated  by  using  Eqs.  (11)  and  (16),  respectively,  and  these  numbers  were  then 
used  in  Eq.  (10).  In  each  case,  the  image  width  cr^  was  taken  to  be  the  short-exposure 
beam  size  value  listed  in  Tables  1-7,  and  the  image  field  of  view  angle  ^  was  also  taken 
to  be  the  short-exposure  beam  size. 
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For  each  scenario,  the  optimum  laser  wavelength  was  chosen  and  then  the  input 
parameters  for  the  lidar  equation  were  chosen  based  on  specifications  of  commercial 
laser,  choices  of  receiver  and  operational  parameters,  and  modeled  or  published  values  of 
atmospheric  characteristics.  The  results  of  these  calculations  are  shown  in  Tables  8-11. 


Table  8.  ABL-ACT  -  50  km  Range 


c  ^ 

RMS  Position  Error  (prad) 

RMS  Diff.  Motion  (prad) 

(day) 

(night) 

lE-17 

2.33E-02 

2.33E-02 

1.1 

lE-16 

8.47E-02 

8.46E-02 

3.5 

lE-15 

4.23E-01 

4.15E-01 

11 

Table  9.  ASTRONOMY 


c  ^ 

RMS  Position  Error  (prad) 

RMS  Diff.  Motion  (prad) 

(night) 

2.20E-13 

4.09E-03 

0.7 

7.00E-13 

9.35E-03 

1.3 

2.20E-12 

1.93E-02 

2.3 

lE-11 

0.05 

5 

Table  10.  ABL  -  100  km  Range 


c  ^ 

RMS  Position  Error  (prad) 

RMS  Diff  Motion  (prad) 

(day) 

(night) 

l.OOE-18 

5.54E-03 

5.54E-03 

0.4 

5.00E-18 

1.32E-02 

1.32E-02 

0.92 

l.OOE-17 

2.08E-02 

2.08E-02 

1.3 

5E-17 

2.9 
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Table  1 1 .  ABL  -  300  km  Range 


c  ^ 

RMS  Position  Error  (prad) 

RMS  Diff.  Motion  (prad) 

(day) 

(night) 

l.OOE-17 

1.67E-02 

1.67E-02 

0.7 

l.OOE-16 

4.57E-02 

4.53E-02 

1.6 

l.OOE-15 

5.49E-02 

5.49E-02 

2.2 

5.00E-15 

0.22 

0.20 

5 

In  Tables  8-11,  the  RMS  position  errors  are  at  most  a  few  percent  of  the  expected 
differential  motions  due  to  turbulence,  during  either  day  or  night.  This  means  that  the 
position  errors  due  to  photon  statistics  will  have  a  small  effect  on  the  measurements  of 
cl .  The  conclusion  of  the  performance  analysis  is  therefore  that  the  proposed  DIM  lidar 

technique  is  practical  in  several  scenarios,  including  the  ABL  Advanced  Concepts 
Testbed,  astronomy,  and  as  a  Tactical  Decision  Aid  (TDA)  for  the  ABL  aircraft. 

3.3  Hardware  and  Software  Development 
3.3.1  3.3.1  Hardware  development 

The  purpose  of  the  field  experiment  was  to  validate  the  DIM  lidar  concept  by  measuring 
the  strength  of  atmospheric  turbulence  with  this  method  and  with  an  independent  sensor. 
The  specific  objectives  were  two-fold;  1)  to  show  that  the  method  accurately  measures 
cl  values  throughout  a  range  of  turbulent  conditions,  and  2)  to  show  that  the  two-way 
(lidar)  technique  gives  the  same  results  as  the  conventional,  one-way  DIM  method. 

To  simplify  the  instrumentation,  we  built  a  hard-target  analog  of  a  DIM  lidar.  A  diagram 
of  this  apparatus  is  shown  in  Fig.3.  The  beam  fi-om  a  25  mW  He-Ne  laser  was  expanded 
with  a  high-quality  optical  system  10  cm  in  diameter  and  focused  on  a  diffuse  target.  The 
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receiver  system  was  based  on  a  40-cm  diameter  Newtonian  telescope.  An  aperture  mask 
with  two  10  cm  diameter  sub-apertures  separated  by  30  cm  was  mounted  in  front  of  the 
telescope  to  select  portions  of  the  incoming  wavefront.  Two  small  glass  plates  were  used 
to  displace  the  images  in  the  image  plane. 

A  lens  assembly  behind  the  focal  plane  formed  a  magnified  image  in  a  CCD  camera  and 
allowed  the  differential  motion  to  be  measured  accurately.  An  interference  filter  between 
the  lenses  rejected  all  but  a  narrow  range  of  wavelengths  centered  on  the  laser  line  and 
improved  the  signal-to-noise  ratio  by  blocking  sunlight.  The  CCD  had  a  650  x  693  pixel 
array  and  a  digital  interface  with  a  Pentium-based  computer  through  the  PCI  bus.  A 
sampling  rate  of  100  Hz  was  achieved  by  reading  a  128  x  64  pixel  sub-array  that  covered 
both  images.  Software  was  developed  to  calculate  the  centroid  position  for  each  image, 
to  compute  the  variance  crj  ,  and  finally  to  calculate  .  For  each  data  point,  1000  data 
frames  were  acquired  to  estimate  the  DIM  variance. 


Interference 

Filter 


Fig.  3.  Diagram  of  the  hard-target  analog  of  the  DIM  lidar  used  in  the  field  test. 
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3.3.2  Software  development 


The  purpose  of  the  specially  developed  software  was  to  acquire  data  and  compute  the 
structure  characteristic  .  The  software  package  included  the  data  acquisition  software, 
the  algorithm  for  determining  the  energy  centroid,  the  algorithm  for  calculating  rms 
differential  motion,  cr^ ,  and  the  algorithm  for  . 

The  frame  rate  was  a  primary  concern  in  the  data  acquisition  software  development.  The 
camera  was  shipped  with  a  simple  software  development  kit  (SDK)  and  C  example 
programs.  The  full  focal  plane  array  of  the  EDC  2000N  camera  consisted  of  650x693 
pixels.  A  pixel  read  time  of  about  2  million  pixels/sec  limited  the  frame  rate  to  about  12 
frames/sec.  In  order  to  read  frames  at  rates  of  80  Hz  -lOOHz,  which  were  necessary  for 
the  plarmed  experiment,  the  sub-array  option  was  employed.  The  sub-array  size  was 
computed  to  encompass  the  expected  separation  of  the  spots  in  the  image  plane  as  well  as 
the  expected  size  of  the  spots  and  their  movement  as  presented  in  Table  1. 

Based  on  the  estimates  of  spot  separations,  spot  size  and  spot  motion  characteristics,  a 
128x64  pixel  sub-array  was  selected.  The  read-out  time  for  this  sub-array  was  about  4 
ms.  A  variable  integration  time  was  selected,  in  the  range  from  1  ms  to  6  ms,  to 
accommodate  variations  of  the  laser  irradiance  and  a  background.  Full  frame  read-out 
time  includes  an  integration  time,  a  mean  overhead  time  (which  ranges  from  2-3  ms  per 
frame),  and  the  pixel  transfer  time.  Depending  on  the  integration  time,  a  frame  rate  in  the 
range  from  80  frames/sec  to  120  frames/sec  was  demonstrated  with  the  sub-array  option. 

The  software  operates  as  follows;  first,  the  imager  is  blocked  and  a  background  of  the  full 
focal  plane  array  is  recorded.  Then,  the  software  captures  a  number  of  frames,  subtracts 
the  background  from  the  frames,  and  computes  the  mean  X  and  Y  values  of  the  centroids 
for  all  of  the  frames.  The  mean  centroid  location  is  used  as  the  center  of  the  analysis 
region.  The  SDK  requires  that  the  sub-region  be  attached  to  the  left  side  of  the  entire 
array.  The  software  can  place  the  analysis  region  anywhere  in  the  image  inside  the  focal 
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plane  array.  However,  the  performance  degrades  as  the  analysis  region  moves  away  from 
the  left  side  of  the  image  because  the  camera  is  transferring  additional  pixels. 


Once  the  alignment  of  the  analysis  region  is  completed,  the  software  collects  a  set  of 
1000  data  frames,  transfers  them  to  memory,  and  subtracts  the  background  from  each 
frame.  The  X  and  Y  centroids  are  then  computed  for  all  of  the  frames.  The  centroid  is 
computed  as  X  =  ,  where  />,.is  the  image  brightness  at  the  i'*’  pixel,  and 

<■  i 

Ij  is  the  distance  from  the  pixel  to  the  center  of  the  sub-array.  From  the  centroids  the 
differences  are  computed:  aij)  =  X^-X,  and  fiU)  =  7,  - 1^ ,  where  j  is  the  number  of 
the  frame.  Once  all  the  differences  have  been  computed,  the  software  computes  the  mean 
and  variance  of  the  differential  image  motion  and  the  structure  characteristic  from  Eq. 
(4). 


3.4  Field  test 


In  order  to  show  that  the  DIM  lidar  technique  gives  correct  values  of  ,  we  tested  the 

hard-target  analog  against  a  conventional  scintillometer  manufactured  by  Scientific 
Technology,  Inc,  which  measures  scintillations  of  an  infrared  light  and  provides  path- 
averaged  readings  of  .  The  test  was  performed  on  a  300  m  horizontal  path  in  an  open 

area  covered  with  a  grass.  The  instrumentation  included  the  hard-target  analog  of  the 
DIM  lidar,  the  scintillometer,  and  weather  instruments.  Weather  instruments  included  a 
wet  and  dry  bulb  hygrometer,  an  ultrasonic  anemometer,  and  an  infrared  temperature 
sensor.  The  hygrometer  provided  air  temperature  and  relative  humidity,  the  anemometer 
measured  wind  velocity  and  direction,  and  the  infrared  sensor  provided  the  ground 
temperature.  To  show  that  the  DIM  variances  on  two-way  and  one-way  paths  coincide 
with  each  other,  we  also  collected  data  by  using  a  conventional  one-way  DIM 
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configuration.  In  this  case,  a  laser  source  at  the  diffuse  target  location  was  used  to 
illuminate  the  DIM  receiver. 

The  DIM  sensor  and  the  scintillometer  were  operated  on  adjacent  paths.  The  diffuse 
target  and  the  scintillometer  transmitter  were  located  at  one  end  of  the  measurement  path, 
and  the  DIM  sensor,  the  scintillometer  receiver,  and  the  data  acquisition  system  were 
located  at  the  other  end  of  the  path.  The  anemometer  was  placed  in  the  middle  range. 
The  data  from  both  optical  instruments  were  collected  at  10-second  intervals.  The  height 
of  the  optical  path  above  the  terrain  was  1.4  m  at  both  ends  and  it  varied  from  0.75  m  to 
2.3  m  at  various  points  along  the  path.  The  average  height  was  \Am. 

The  field  experiment  was  performed  on  October  5  -  7,  1999.  On  October  7,  the  data 
were  continuously  collected  from  0845  to  1245  EDT.  The  sky  was  clear  initially,  with 
some  cmnulus  clouds  toward  the  end  of  the  test.  The  air  temperatiure  increased  from  1 5 
to  22  ®C,  and  the  grass  was  dew-covered  for  most  of  the  test.  The  wind  was  calm  at  first, 
changing  to  a  light  breeze  (<5  m/s). 


3.5  Analysis  of  the  data 

The  data  collected  in  the  test  were  processed  and  analyzed  by  using  the  specially 
developed  software.  Figure  4  shows  a  time  history  of  the  October  7  measurements.  The 
C]  values  were  measured  with  the  DIM  sensor  in  both  two-way  and  one-way 
configurations,  and  also  with  the  scintillometer.  The  results  are  plotted  as  two-minute 
averages  for  clarity.  Figure  4  shows  that  the  Cl  values  measured  with  the  hard-target 

analog  of  the  DIM  lidar  and  with  the  scintillometer  agree  well  with  each  other.  It  also 
shows  that  the  data  collected  with  the  one-way  and  two-way  DIM  configurations  are 
consistent  with  each  other  and  that  they  both  agree  with  the  scintillometer  data. 

A  scatter  plot  of  the  10  second  data  collected  with  the  hard-target  analog  and  the 
scintillometer  is  shown  in  Figure  5.  This  also  shows  good  agreement  between  the 
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readings  from  the  two  sensors.  The  slope  of  the  straight  line  obtained  from  a  linear 
regression  analysis  is  1.2  and  the  regression  coefficient  is  0.77.  The  differences  between 
the  readings  from  the  DIM  sensor  and  scintillometer  could  be  partly  due  to  variations  of 
along  the  path.  The  path-weighting  function  for  the  DIM  sensor  is  maximal  near  the 

DIM  receiver.  At  this  point,  the  height  above  the  ground  is  1.4  m.  The  path-weighting 
function  for  the  scintillometer  is  maximal  at  the  middle  range,  where  the  height  is  about 
2.3  m.  The  well-known”  decrease  of  with  height  would  be  expected  to  lead  to  a 

difference  between  the  readings  from  the  two  sensors.  Different  local  albedos  and 
amounts  of  surface  moisture  at  the  near  end  of  the  path  and  at  the  middle  range  could  also 
contribute  to  the  observed  differences  between  the  values. 


Hour  of  Day  (EDT) 


Fig.  4.  Time  histories  of  the  values  measured  with  a  commercial  scintillometer 
and  the  hard-target  analog  of  the  DIM  lidar.  Curve  1:  scintillometer.  Curve  2: 
hard-target  analog.  Curve  3:  one-way  configuration  of  the  DIM  sensor. 


The  measurements  reported  here  were  performed  throughout  a  range  of  turbulence 
conditions.  The  strength  of  turbulence  is  usually  characterized  by  the  Rytov 
variance, =  1.23A:’^®CfZ"'* .  For  the  wavelength  of  A  =  0.63  10"*/m  ,  distance 
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L  =  300/w ,  and  <3-10  the  Rytov  variance  was  cr^  <  2 ,  so  the  validation  test 

was  performed  in  the  regime  of  weak  and  intermediate  scintillation. 

The  experimental  results  validate  the  concept  of  the  DIM  lidar  and  show  that  the  structure 
characteristic  can  be  measured  accurately  with  this  technique.  The  theoretical  and 

experimental  underpinnings  are  now  in  place  to  develop  a  true  DIM  lidar  that  will 
employ  a  pulsed  laser,  atmospheric  backscatter,  and  a  range-gated  receiver  to  measure 
profiles  of  cl . 


log  from  scintillometer 


Fig.  5.  Scatter  plot  of  the  measurements  with  the  hard-target  analog  of  the  DIM 
lidar  versus  the  scintillometer  data. 
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3.6  Scientific  paper 


A  scientific  paper  was  written  and  submitted  to  Optics  Letters.  The  paper  was  accepted 
for  publication. 

3.7  Effects  of  laser  speckles  and  scintillations 

Air  Force  Experience  with  an  astronomical  Differential  Image  Motion  Monitor  (DIMM) 
at  the  Starfire  Optical  Range  (SOR)^^  and  at  North  Oscura  Peak  (NOP)^^  has  shown  that  a 
speckle  phenomenon  degrades  the  performance  of  the  sensor  during  conditions  of  strong 
turbulence.  The  primary  purpose  of  the  analysis  presented  in  this  section,  which  was 
documented  in  a  Technical  Memorandum  and  sent  to  the  sponsor,  was  to  assess  the  effect 
of  this  phenomenon  on  the  performance  of  our  proposed  DIM  lidar. 

3.7.1  Laser  speckle  phenomenon 

It  is  known  that  for  some  turbulent  conditions  a  star  image  breaks  into  multiple  speckles. 
This  phenomenon  occurs'^’’^  when  the  telescope  diameter,  D  ,  exceeds  the  Fried 
parameter,  ro ,  D/kq  >1.  The  speckles  can  introduce  a  measurement  error  in  the  centroid 
position  for  both  open-loop  tilt  measurements  for  astronomical  observations  and  for 
sensing  of  ro  with  a  DIM  sensor.  A  simulation  performed  by  a  researcher  at  SOR 
revealed  that  the  maximum  error  in  the  centroid  position  is  introduced  when  the  number 
of  speckles  in  the  image  plane,  »  (D  /  ,  is  severely  limited,  or  D  /  «  1  -  3 .  The 

error  decreases  with  an  increasing  number  of  speckles,  N,  because  multiple  speckles 
compensate  each  other.  According  to  the  simulation  performed  by  Capt.  J.  Barchers,  the 
speckle-induced  relative  error  in  the  centroid  position  for  D  /  Tq  »  1  -  3  is  in  the  range  of 

15-20  %.  However,  this  estimate  is  not  applicable  to  the  DIM  lidar  because  the 
phenomenon  observed  at  SOR  and  NOP^^  only  occurs  with  point  sources  such  as  a  star  or 
a  distant  laser  and  does  not  occur  with  the  DIM  lidar  since  the  scattering  volume 
illuminated  by  a  laser  beam  is  not  a  point  source. 
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3.7.2  Coherent  point  source 


In  the  conventional  DIM  technique,  a  receiver  is  illuminated  by  a  coherent  point  source 
(a  star  or  a  laser).  If  the  initial  radius  of  the  beam,  a^,  satisfies  the  condition 

Q  =  kal  / L«l,  then  the  source  is  a  point  source.  Here  k  =  2n ! X  is  the  wavenumber 

and  L  is  the  measurement  range.  The  spatial  coherence  of  an  optical  wave  is  degraded 
exclusively  by  turbulence.  For  a  plane  wave  the  mutual  coherence  function  is  expressed 
through  the  corresponding  phase  structure  function,  , 

=  exp(-^D^(p))  =  exp(-(yo//?o)''') ,  (17) 


where  is  the  coherence  radius,  =  (1  ASk^LCl  is  the  refi-active  index 

structure  characteristic.  This  scale  relates  to  the  Fried  parameter  by  rp  =  2.1yOo .  For  a 
spherical  wave,  or  point  source,  the  coherence  radius  equals 

L 

Pc={\A5k^  (18) 

0 

and  Pc  =  for  cl  =  const. 

In  the  case  of  a  DIM  lidar,  a  transmitted  beam  illuminates  a  scattering  volume  which 
contains  a  large  number  of  small  particles,  or  molecules,  randomly  distributed  throughout 
the  volume.  The  separation  of  the  particles  is  generally  much  greater  than  an  optical 
wavelength.  As  a  consequence,  the  scattered  radiation  from  one  scatterer  will  be 
completely  incoherent  with  respect  to  that  from  any  other.  This  mechanism  is 
identical  to  that  for  a  diffuse  target,  or  extended  incoherent  source.  The  longitudinal  size 
of  the  scattering  volume  is  determined  by  the  range  gate  of  the  imager,  which  varies  from 
several  himdreds  of  meters  to  several  kilometers.  The  transverse  size  of  the  volume 
coincides  with  the  short-exposme  size  of  the  beam,  and  typical  values  of  the  Fresnel 
parameter  for  the  scattering  volume  size  are  typically  greater  than  unity,  Q  >  1. 
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Due  to  random  phase  perturbations  at  different  positions  of  the  scatterers  under  the 
reflection  of  the  coherent  light  from  the  multiple  scatterers  a  speckle  pattern  in  the  image 
plane  occurs, but  this  is  a  different  phenomenon  from  the  stellar  image  case,  and  the 
speckle  scales  are  much  smaller.  This  is  shown  in  the  next  section 


3.7.3  Extended  incoherent  source 

References  16  and  17  give  the  mutual  coherence  function  of  the  light  scattered  by 
multiple  scatterers  in  the  turbulent  atmosphere  as 

where  is  the  mutual  coherence  function  of  a  laser  beam  scattered  in  a  homogenous 
medium  (C^  =0),  and  is  the  coherence  function  of  a  spherical  wave  in  a  turbulent 
atmosphere.  The  function  Pi''?  has  the  form 


■p(a) 

^(2) 


exp[-^  +  ^R^], 
Pea  L 


(20) 


where  =Llka^  is  the  coherence  radius  determined  by  the  van  Cittert-Zemike 
theorem,  the  effective  radius  of  the  transmitted  beam,  p  is  the  backscatter 

coefficient,  and  Iq  is  the  laser  intensity.  In  most  cases  the  coherence  radius,  p^ ,  is  much 
smaller  than  the  turbulent  coherence  scale,  «  p^.  As  a  consequence,  the  first  term  in 
equation  (12)  is  predominant,  and  the  coherence  scale  of  the  scattered  light  equals  p^^ . 
Consequently,  the  size  of  the  image,  which  is  determined  by  diffraction  of  the  coherence 
scale.  Pea,  ai^f  !  kp^a ,  where/is  the  focal  length,  coincides  with  the  angular  size  of  the 

transmitted  beam,  =aj  f  =  a^lL.  Thus,  the  size  of  the  image  is  determined  by  the 
angular  size  of  the  beam  ,  not  by  turbulence  as  for  a  point  source. 
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Now  we  estimate  the  number  of  speckles  in  the  image  spot  of  a  DIM  lidar  for  various 
scenarios.  In  the  proof-of  concept  test  in  Phase  I,  a  CW  laser  with  X  =0.63  fjm  was  used. 
The  initial  diameter  of  the  Gaussian  beam  was  2a^  =  6.4  cm,  and  the  receiver  diameter 

was  Ior  =10  cm.  The  diffractive  size  of  the  beam  focused  at  a  distant  diffuse  target, 
=1/  ka^,  is  <pd=  3.1  fjrad.  However,  due  to  a  multi-mode  structure  of  the  laser,  the 

actual  beam  size  was  about =  4x^^.  Since  for  the  aperture  of  Iur  =  10  cw  the 
speckle  size,  ,  is  ^5=  1.8  fjm,  the  number  of  speckles  in  the  reflected  spot  is 

^  =  The  number  of  speckles,  N,  increases  with  the  strength  of 

turbulence  because  turbulent  broadening  of  the  beam  increases  (p^ . 

Figures  6  a-c  show  images  collected  in  preliminary  tests  of  the  hard-target  analog  of  a 
DIM  lidar.  The  images  were  collected  in  the  hallway  over  a  50  m  path  without  a  heater 
(6  a),  this  case  corresponds  to  very  weak  turbulent  conditions,  over  a  50  m  path  with  a 
heater  (6  c),  this  case  corresponds  to  very  strong  turbulence,  and  over  150  m  outdoors 
during  daytime  (6  b).  The  last  case  corresponds  to  moderate  turbulent  conditions.  It  is 
seen  that  for  all  conditions  the  number  of  speckles  is  large,  A^>1,  therefore  the  error  m 
the  centroid  position  caused  by  speckles  is  expected  to  be  small.  Numerical  simulation  of 
a  single  shot  image  for  a  coherent  lidar  (see  Ref.  8,  Fig  2)  provides  a  speckle  pattern 
which  is  similar  to  those  shown  in  Figures  6  a  and  b. 


Figure  6  a-c.  Images  collected  with  a  hard-target  analog  of  the  DIM  lidar  in  very 
weak  (a),  moderate  (b),  and  very  strong  turbulence  condition  (c). 

Considering  an  eye  safe  uv  lidar,  the  maximum  measurement  range  at  the  horizontal  path 
at  NOP  is  limited  by  the  transmittance  of  the  UV  light  to  about  L  =20  km.  For  L  =20  km 
and  a  maximum  value  of  Cl  =2.1 m-^'\  the  short-exposure  beam  size  at  NOP  is 
(Peff  ,  therefore  the  number  of  speckles  is  V=118  for  2a«=  16  cm,  and  it  is 

N—2'i6  for  2aR=  32  cm.  For  the  astronomical  scenario  with  a  maximum  range  of  L=  20 
km  and  zenith  angle  of  45°,  N=  95  for  2a^  =  16  cm,  and  V=190  for  Iur  =  32  cm. 


The  long-exposure  beam  radius  was  estimated  from  Eq.(7),  the  estimates  for  the  short- 
exposure  size  of  the  beam  were  obtained  from  Eq.(8),  and  the  one-axis  beam  wander 
variance  was  computed  from  Eq.  (9).  From  equation  (7)  it  is  seen  that,  for  long  paths  and 
strong  turbulence,  when  the  contribution  of  the  turbulent  term  is  predominant, 
^eff  **  (^n )  ’  the  nxunber  of  speckles  in  the  reflected  spot  increases  with  the  strength  of 
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turbulence  as  iV  «  (C^  thus  reducing  the  speckle-induced  error  for  the  centroid.  In  the 
next  section  we  will  quantify  this  effect 


3.7.4  Error  analysis 

To  quantify  the  measurement  error  in  the  centroid  position  caused  by  laser  speckles,  the 
centroid  variance  in  a  speckle  image  is  calculated  in  the  absence  of  turbulence.  A 

centroid  is  given  by  C  =  P~'  j^RI(R)d^R ,  where  P  is  the  flux  through  the  aperture,  /  (R) 

is  the  intensity  distribution  in  the  image  plane,  and  S  =  ml  is  the  area  of  the  receiving 
aperture.  By  using  this  expression,  the  centroid  variance  cr^  can  be  expressed  through  the 
intensity  covariance  given  by 

al  =  P-^  lld^Rd^pA(R)Aip)(R^  - p^  /4)(/(^  +  plT)I{R  -pH)).  (21) 

Here  A  (R)  is  the  aperture  function  of  the  receiver,  and  the  angular  brackets  denote 
ensemble  averaging. 

Due  to  random  phase  perturbations  of  the  scattered  light,  the  optical  field  fi-om  a 
scattering  volume  or  a  diffuse  target,  is  a  complex  Gaussian  random  variable. The 
fourth  moment  of  the  random  Gaussian  variable  relates  to  the  second  moment  by  the 
equation,^® 


(7, /,>  =  |r,p +(/>“,  (22) 

where  r(2)  is  given  by  equation  (12).  By  using  equations  (19)  and  (22),  as  well  as  a 
Gaussian  approximation  for  the  aperture  transmission  function  (y4(i?)  =  exp(-i?^  Hi)), 
one  can  obtain  g\=  p],  where  Pj  is  the  intensity  correlation  scale.  From  equation  (22) 
it  follows  that  the  intensity  correlation  scale  is  given  by  pj  =  p^^  /  . 
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Since  the  structure  eharaeteristic,  ,  is  expressed  through  the  mean  square  difference 
((«,  where  aj  and  «2  are  the  centroids  of  the  two  images,  the  absolute  error  for 

Cl  caused  by  the  speckle  phenomenon  equals  double  the  error  for  the  centroid,  or 
ACl  =2<jI  .  We  estimate  the  relative  error,  e  =  ACl  /  Cl ,  for  the  case  L  >  lOkm ,  when 
the  effective  size  of  the  beam  is  determined  by  the  turbulence,  or  the  third  term  in 
equation  (7).  The  estimates  show  that  the  rms  differential  motion,  ,  in  this  regime 
relates  to  the  short-exposure  beam  radius  by  /y ,  where  yis  in  the  range  from 

1  to  1.7.  Since  Cl  ~  ,  the  relative  error,  s  =  ACl  /Cl ,  equals 

£  =  2al/alj^  =r^ /N .  (23) 

As  the  number  of  speekles  for  2aR=  16  em  is  N=1 1 8,  the  relative  error  for  is  ^  ^  3% 
For  the  receiver  with  2aR=  32  cm  the  number  of  speckles  is  N=236,  and  the  relative  error 
is  £  <1.5%  The  implication  is  that,  due  to  the  large  number  of  speckles  in  the  reflected 

(~i2 

spot,  the  measurement  error  for  the  eentroid  position  and  consequently  for  measured 
with  the  DIM  lidar  is  small. 

3.7.5  Laser  scintillations 

It  is  well  known  that  the  planets  twinkle  weaker  than  stars  situated  at  the  same  zenith 
angle.^’  Tatarskii**  studied  this  phenomenon  theoretieally  and  showed  that  if  the  angular 
size  of  the  incoherent  source  exceeds  the  angular  size  of  the  Fresnel  zone,  >  -yjA/L  , 
then  the  scintillations  are  averaged  over  the  size  of  the  source.  Aceording  to  Tatarskii’s 
estimates,*'  the  scintillation  variance  reduces  by  a  factor  of  about  3  when 
(pg^  l-JXI L  >1. 

Researchers  from  NOAA  used  the  same  phenomenon  to  develop  a  saturation-resistant 
optical  scintillometer.^*  In  this  scintillometer,  both  the  incoherent  source  and  the  receiver 
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exceed  the  radius  of  a  Fresnel  zone  for  the  optical  path,  VaT  .  It  was  shown  that,  for 
such  a  source  and  receiver,  there  is  no  saturation  effect  even  in  the  strong  turbulence 
regime,  and  the  calibration  remained  unchanged  in  both  the  weak  and  strong  scintillation 
regime.  Experimental  validation  performed  over  an  L=500m  path  in  strong  turbulence, 
when  the  Rytov  variance,  =  =  25,  confirmed  this  conclusion. 

As  the  size  of  the  scattering  volume  and  the  receiver  in  the  DIM  lidar  exceeds  the  Fresnel 
zone,  a  similar  saturation  resistance  capability  is  expected.  Indeed,  for  the  moderate 
strength  of  turbulence  at  NOP,  Cf  which  corresponds  to  =10c/n  at 

X  =  1.06//m ,  the  R5dov  variance  is  Pi  =  1.6.  However,  the  effective  size  of  the  beam  for 
X  =  0.2>5\fMn  and  D=  20  km  exceeds  the  angular  size  of  the  Fresnel  zone, 
■\IX/ L  -  Ajurad ,  and  the  receiver  with  an  aperture  diameter  of  2ajj=  32  cm  is  twice  as 
large  as  the  Fresnel  zone.  For  this  reason  one  can  expect  that,  due  to  averaging  of 
intensity  fluctuations  over  the  scattering  volume  and  the  receiver,  the  scintillations  will 
be  reduced  to  the  level  of  weak  turbulence. 

For  the  maximvun  strength  of  turbulence  at  NOP,  C^=2.7-10~'®m’^^\  which 
corresponds  to  Tq  =  4c7m  at  A  =  1 .06//»i ,  the  Rytov  variance  over  an  Z  =  20  km  path  is 
Pi  =  7.  However,  the  effective  size  of  the  beam,  ,  increases  with  the  strength  of 

turbulence.  As  a  consequence,  for  these  conditions,  the  size  of  the  beam  is 
» 1.6ijrad ,  and  it  exceeds  the  Fresnel  zone  by  a  factor  of  2.  Due  to  averaging  of  the 

scintillations  over  the  size  of  the  scattering  volume  and  the  receiver,  according  to  the 
preliminary  estimates,  one  can  expect  that  no  saturation  effect  will  be  observed.  The 
implication  is  that  a  DIM  lidar,  similar  to  the  NOAA  scintillometer,  would  be  expected  to 
be  saturation  resistant. 

These  estimates  also  show  that  for  inclined  propagation  paths  for  astronomical  scenarios, 
the  weak  scintillation  regime  will  be  observed  during  either  day  or  night  even  for  the 
strongest  turbulence  near  the  groimd,  if  the  zenith  angle  of  the  measurement  path  is 

e<6Q\ 
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3.7.6  Sky  background 


The  high  sky  background  level  is  responsible  in  part  for  the  failure  of  a  conventional 
DIM  sensor  to  measure  ro  during  daytime  at  SOR.  In  the  DIM  lidar  this  difficulty  can  be 
overcome  through  the  use  of  a  narrow  band  filter  centered  at  the  laser  wavelength.  This 
will  provide  a  high  signal-to-  noise  ratio  for  both  day  and  night  operations. 

3.7.7  Design  consideration 

Another  possible  reason  for  the  failure  of  the  conventional  DIM  sensor  in  the  presence  of 
strong  turbulence  is  that  the  size  of  the  spots  in  the  image  plane  and  their  motion  could 
exceed  the  FOV,  and/or  two  images  could  overlap.  An  accurate  design  of  the  DIM  lidar 
using  the  flow-down  design  approach  will  allow  for  an  accommodation  of  the  expected 
spot  size  and  motion  for  all  operational  conditions  and  propagation  scenarios.  The  proof- 
®f"Concept  test  with  a  hard-target  analog  of  the  DIM  lidar  confirmed  that  even  for  the 
strongest  turbulence  the  images  of  two  spots  are  within  the  FOV  and  do  not  overlap. 

3.7.8  Summary 

The  analysis  performed  leads  to  the  following  conclusions: 

1.  The  phenomenon  observed  at  SOR  and  NOP,  which  causes  the  images  to  break  up 
into  individual  speckles,  only  occius  with  point  sources  such  as  a  star  or  a  distant 
laser.  Because  the  scattering  volume  of  the  DEM  lidar  will  be  an  extended  source, 
this  particular  phenomenon  will  not  occur,  regardless  of  the  strength  of  turbulence. 

2.  Using  a  laser  as  a  source  of  light  will,  however,  cause  a  different  kind  of  speckle 
pattern  to  appear  in  the  images,  due  to  the  scattering  of  coherent  light  by  multiple 
scatterers.  This  is  different  fi'om  the  point  source  phenomenon,  and  the  scale  of  these 
speckles  is  much  smaller.  For  the  scenarios  and  system  parameters  of  interest  here, 
there  will  typically  be  on  the  order  of  100  or  more  speckles  in  each  image.  This  will 
cause  a  r.m.s.  measurement  error  in  the  centroid  positions  on  the  order  of  1  percent. 
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regardless  of  turbulence  conditions.  Examples  of  images  from  our  field 
demonstration  apparatus  that  show  this  small-scale  speckle  pattern  are  included  in 
this  report. 

3.  Strong  turbulence  will  cause  scintillations  in  the  laser  beam,  and  this  could,  in 
principle,  make  measurements  difficult  due  to  large  fluctuations  in  signal  levels,  but 
the  sizes  of  the  scattering  voliune  and  receiver  aperture  for  the  DIM  lidar  will  tend  to 
average  out  such  effects. 

4.  The  DIM  lidar  will  employ  a  narrow-band  optical  filter  to  eliminate  most  of  the  sky 
background  and  maintain  a  good  signal-to-noise  ratio  during  daytime  operation.  This 
is  possible  because  the  DIM  lidar  will  use  a  laser  as  a  source  of  illumination,  and  a 
laser  is  a  narrow-band  source.  Thus,  the  use  of  a  narrow-band  optical  filter  will 
eliminate  most  of  the  sky  background  and  provide  a  good  signal-to-noise  ratio  during 
daytime  operation. 

4.0  CONCLUSIONS 

All  the  tasks  and  milestones  of  Phase  I  were  accomplished.  The  following  results  were 

obtained: 

1.  The  conceptual  design  and  performance  analysis  of  the  DIM  lidar  were  performed 
and  the  theoretical  imderpinnings  of  the  DIM  lidar  were  developed.  It  was  shown  that 
the  differential  image  motion  of  a  lidar  return  for  the  propagation  scenarios  of  interest 
can  be  accurately  measured; 

2.  The  hardware  of  the  hard-target  analog  of  the  DIM  lidar  was  designed  and  built,  and 
the  software  was  developed; 

3.  A  proof-of-concept  test  was  performed.  The  experimental  results  validate  the  DIM 
lidar  concept  and  show  that  the  structure  characteristic  can  be  accurately  measured 
with  this  technique; 
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4.  The  effects  of  the  laser  speckles  and  scintillations  on  the  performance  of  the  DIM 
lidar  were  analyzed.  It  was  shown  that  the  error  in  the  centroid  position  caused  by  the 
laser  speckles  does  not  exceed  a  few  percent.  The  DIM  lidar  has  a  saturation  resistant 
capability,  which  is  due  to  averaging  of  the  scintillations  over  the  scattering  volume 
and  the  receiving  aperture.  The  use  of  a  narrow-band  filter  will  eliminate  most  of  the 
sky  background  and  provide  a  good  signal-to-noise  ratio  during  daytime  operation; 

5.  A  scientific  paper  was  written  and  submitted  for  publication  to  Optics  Letters. 
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